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Abstract 

A general analysis of the sensitivities of neutron beta-decay experiments to manifestations of 
possible interaction beyond the Standard Model is carried out. In a consistent fashion, we take 
into account all known radiative and recoil corrections arising in the Standard Model. This provides 
a description of angular correlations in neutron decay in terms of one parameter, which is accurate 
to the level of ~ 10~^. Based on this general expression, we present an analysis of the sensitivities to 
new physics for selected neutron decay experiments. We emphasize that the usual parametrization 
of experiments in terms of the tree level coefficients a, A and B is inadequate when the experimental 
sensitivities are at the same or higher level relative to the size of the corrections to the tree level 
description. 
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I. INTRODUCTION 



The relative simplicity of the decay of the free neutron makes it an attractive laboratory 
for the study of possible extensions to the Standard Model. As is well known, measurements 
of the neutron lifetime and neutron decay correlations can be used to determine the weak 
vector coupling constant, which, in turn, can be combined with information on strange par- 
ticle decay to test such notions as the universality of the weak interaction or to search for (or 
put a o„) „o„.a.da.d coup,,,,. ,ee. to. e.au,ple. QBflQQflBa aud .efe.euce. 
therein). It is less widely appreciated that precision measurements of the correlations in 
neutron decay can, in principle, be used as a test of the standard model without appeal to 
measurements in other systems. In particular, the detailed shape of the decay spectra and 
the energy dependence of the decay correlation are sensitive to non-standard couplings. The 
extraction of such information in a consistent fashion requires a rather delicate analysis, as 
the lowest order description of the correlation coefficients (and their energy dependencies) 
must be modified by a number of higher order corrections that are incorporated within the 
Standard Model. These include such effects as weak magnetism and radiative corrections. 
Rece^ly Q effective «e.d theory Ka. bee,, used to „,co,.po,.ate all .ta„da,.d ,„odel effects ,n 
a consistent fashion in terms of one parameter with an estimated theoretical accuracy on the 
order of 10^^. Because this accuracy is well below that anticipated in the next generation of 
neutron decay experiments (see, for example, papers in jlo|), this analysis provides a useful 
framework for the exploration of the sensitivity of various experiments to new physics. 

In this paper, we extend the description of neutron beta-decay of by including the 
most general non-standard beta-decay interactions. Our framework provides a consistent 
description of the modifications of the beta-decay observables at a level well below that 
anticipated in the next generation of experiments. Not surprisingly, we find that the different 
experimental observables have quite different sensitivities to the form of hypothetical non- 
standard couplings (i.e. vector, scalar, etc.). 
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II. NEUTRON /3-DECAY BEYOND THE STANDARD MODEL. 



The most general description of neutron /3-decay can be done in terms of low energy 



constants Cj by the Hamiltonian 



11 



Hint = {'lpp'lpn)iCs'lpe'lpu + C's'lpe'J^i'u) 

- {'4^pltilb'4^n){CAi^eli,lb'4^u + C'^^'eT^^'i^) 

+ (V'p75V'n)(Cp^e75V'i^ + C'p^e^i^) (l) 

+ Hermitian conjugate, 

where the index i = V , A, S^T and P corresponds to vector, axial-vector, scalar, tensor and 
pseudoscalar nucleon interactions. In this presentation, the constants Ci can be considered 
as effective constants of nucleon interactions with defined Lorentz structure, assuming that 
all high energy degrees of freedom (for the Standard model and any given extension of the 
Standard model) are integrated out. In this paper we consider only time reversal conserving 
interactions, therefore the constants Ci can be chosen to be real. (The violation of time 
reversal invariance in neutron decay at the level of considered accuracy would be a clear 
manifestation of new physics and thus does not require an analysis of the form contained 
here.) Ignoring electron and proton polarizations, the given effective Hamiltonian will result 



in the neutron /3-decay rate [l^ in the tree approximation (neglecting recoil corrections and 
radiative corrections) in terms of the angular correlations coefficients a. A, and B\ 



x(l + b— + a ^ +A — — + B — (2) 

Here, a is the neutron spin; me is the electron mass. Eg, E^, pe, and p,y are the energies 
and momenta of the electron and antineutrino, respectively; and Gp is the Fermi constant 
of the weak interaction (obtained from the /i-decay rate). The function $(-Ee) includes 
normalization constants, phase-space factors, and standard Coulomb corrections. For the 
Standard model the angular coefficients depend only on one parameter A = —Ca/Cv > 0, 
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the ratio of axial- vector to vector nucleon coupling constant (Cy = Cy and Ca = C'j^: 

^--2^, B^2^_. (3, 



1 + 3A2' 1 + 3A2' 1 + 3A2' 

(The parameter h is equal to zero for vector - axial-vector weak interactions.) 

As was shown in the existence of additional interactions modifies the above expressions 
and can lead to a non-zero value for the coefficient h. To explicitly see the influence of a 
non-standard interaction on the angular coefficients and on the decay rate of neutron one 
can re-write the coupling constants Ci as a sum of a contribution from the standard model 
Cf*^ and a possible contribution from new physics dCi. 
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(4) 



We neglect the pseudoscalar coupling constants since we treat nucleons nonrelativisti- 
cally. Defining the term proportional to the total decay rate in eq.® as ^ = (1 + 3A^) one 
can obtain corrections to parameters S,, a, 6, A and B due to new physics as Sa, 6b, 6A 
and 6B, correspondingly. Then, using results of 

5^ = C^^'{6Cv + SC{,) + {5Cv^ + SC{,^ + 6Cs^ + 6C's^)/2 

+ 3[C^^{6Ca + SC'a) + i6CA^ + 5C'/ + 5Ct^ + 5C't)/2], 
i5h = VT^^[C^^\6Cs + SC's) + 6CsSCv + 6C's6C{, 

+ 3(C^^'(5Ct + 5C^) + 5Ct5Ca + SC'tSC'^)], 
^5a = C^^'i5Cv + 5C'y) + (dCy^ + 5C'y^ - 5Cs^ - SC's^)/2 

- C^^\SCa + 5C'a) - (SCa^ + SC'/ - 5Ct^ - 5C^')/2, 
CSA = -2C^^\6Ca + SC'^) + SC'JC'^ - bCybC'^ 

- [C^^((5C^ + bC'^) + Cf \bCy + bC'y) + bCybC'^ + bC'ybCA - SCsSOr - ^C's^Ct\ 



e 



— 5Cs5C'j^ — dCxSCg — 6Cv6C'j^ — dCA^C'y. (5) 

It should be noted that we have neglected radiative corrections and recoil effects for the 
new physics contributions, because these are expected to be very small. However, Coulomb 
corrections for the new physics contributions are taken into account since they are important 
for a low energy part of the electron spectrum. 

From the above equations one can see that contributions from possible new physics to 
the neutron decay distribution function is rather complicated. To be able to separate new 
physics from different corrections to the expression (j2I), obtained in the tree level of approx- 
imation, one must describe the neutron decay process with accuracy which is better than 
the expected experimental accuracy. Assuming that the accuracy in future neutron decay 
experiments can reach a level of about 10^"^ — 10~^, we wish to describe the neutron decay 
with theoretical accuracy by about 10~^ and our description must include all recoil and 



radiative corrections 



uracy by about iU ana our aescr 

i 13, Q Q, 0, 0, Q, Q 0, Q • 



To do this we will use recent re- 



sults of calculations jsf of radiative corrections for neutron decay in the effective field theory 
(EFT) with some necessary modifications. The results of can be used since they take 
into account both recoil and radiative corrections in the same framework of the EFT with 
estimated theoretical accuracy which is better than 10^^. However, the EFT approach does 
not provide all parameters but rather gives a parametrization in terms of a few (two, in 
the case of neutron decay) low energy constants which must be extracted from independent 
experiments. Therefore, the neutron /3-decay distribution function is parameterized in terms 
of one unknown parameter (the second parameter is effectively absorbed in the axial vector 
coupling constant). If this parameter would be extracted from an independent experiment, 
it gives a model independent description of neutron beta-decay in the standard model with 
accuracy better than 10~^. A rough estimate of this parameter based on a "natural" size 
of strong interaction contribution to radiative corrections gives an accuracy for the expres- 
sions for the rate and the angular correlation coefficients which is better than 10"^ (see j^). 
We vary the magnitude of this parameter in a wide range for the given numerical analysis 
and show that variations of the parameters in the allowed range do not significantly change 



our results at a level well bellow 10~^. Also, unlike [9|, we use the exact Fermi function 
for numerical calculations to take into account all corrections due to interactions with the 
classical electromagnetic field. This gives us the expression for neutron decay distribution 
function as 



dEedilp^dflp^ 



(27r) 



5 \pe\E,{Er''-E,yF{Z,E, 



EeEjj 



Pe-Pu 

^ EeEy 
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where the energy dependent angular correlation coefficients are: 

ME,) = {1 + 3X^) (l + ^^i^^ + ^ e^) 
2 



m 



A(/iy + A)-^ + A(/iy + A)E: 
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(7) 

(8) 
(9) 

(10) 
(11) 



(12) 
(13) 



Here Cy is the finite renormalized low energy constant (LEG) corresponding to the "inner" 
radiative corrections due to the strong interactions in the standard QCD approach; F(Z, Ef,) 
is the standard Fermi function; and the functions 5^'^ and 5^'> are: 
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13 p Vl + /^. 
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where (3 = Pe/Ee- The only unknown parameter Cy is chosen to satisfy the estimate [2^ 
for an "inner" part of the radiative corrections: ^ = 0.02. In Eq.® the custom of 
expanding the nucleon recoil correction of the three-body phase space has been used. These 
recoil corrections are included in the coefficients fi, i = 0,1, ■■■,6 defined in the partial 
decay rate expression, Eq.®. It should be noted that the expression for /2 is an exclusive 



three-body phase space recoil correction, whereas all other /j 



0, 1, 3, ■ ■ ■ , 6 contain a 



mixture of regular recoil and phase space (l/m-Ar) corrections. 

The above expression presents all contributions from the Standard model. Therefore, the 
difference between this theoretical description and an experimental result can only be due 
to effects not accounted for by the Standard model. From the eqs.© we can see that the 
only contributions from new physics in neutron decay are: 



fo{Ee 
fl{Ee 
UEe 
ME, 



fo{E,)+5^- 
fi{E,)+6a, 
ME,)+6A, 
ME,)+6B, 



(16) 



Since possible contributions from models beyond the Standard one are rather complicated, 
we have to use numerical analysis for calculations of experimental sensitivities to new physics. 

III. THE ANALYSIS OF THE EXPERIMENTAL SENSITIVITY TO NEW 
PHYSICS 



To calculate the sensitivity of an experiment with a total number of events N to the 
parameter q we use the standard technique of the minimum variance bound estimator (see. 



for example 



22,|23|). The estimated uncertainties provided by this method correspond to 



TABLE I: Relative statistical error (K) of the standard experiments to different types of inter- 
actions from new physics (Cj constants) provided that these constants have the same values of 
1 • 10-3. 



Interactions 
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A 
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V 


5.26 


3.60 


6.95 


A 


1.73 


1.90 


1.91 


T 


2.59 


7.25 


1.50 


S 


8.70 


26.70 


1.46 


V + A 


2.01 


1.58 


3.86 



one sigma limits for a normal distribution. The statistical error (variance) ag of parameter 
q in the given experiment can be written as 

K 



On 



where 



K- 



Jw{x) 



N 



1 dw(x) 
w{x) dq 



dx 



(17) 



(18) 



/ w{x)dx 

Here w{x) is a distribution function of measurable parameters x. We can calculate the 
sensitivity of the experiment to a particular coefficient Ci or to a function of these coefficients. 
The results for these integrated sensitivities for each type of interaction (Cj) and for the left- 
right model are given in the table for the standard experiments measuring a, A and B 
coefficients in neutron decay, assuming that all coefficients Ci have the same value of 1 ■ 10^^. 
The numerical test shows that results for the coefficients K can be linearly re-scaled for the 
parameters Ci in the range from 10~^ to 10~^ with an accuracy of better than 10%. We 
can see that different experiments have different sensitivities (discovery potentials) for the 
possible manifestations of new physics. 

The given description of neutron /3-decay experiments in terms of low energy constants 
related to the Lorentz structure of weak interactions is general and complete. All mod- 
els beyond the Standard one (new physics) contribute to the Ci values in different ways. 
Therefore, each model can be described by a function of the Ci parameters. To relate these 
C-coefficients explicitly to the possible models beyond the Standard one we can use the 
parametrization of reference [3]. It should be noted that the definitions of reference 12 1 
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used for the Ci coefficients are the same as in 0, except for the opposite sign of Cy, C^, 
and Ca- Therefore, we can re- write the relations of the 5Cj, which contain contributions to 
the Ci from new physics, in terms of the parameters (iji and Aji defined in the paper 0] as: 

5Cy = Cy'^^aii + a^ji + ajiL + ajiji), 

SC'y = —Cy'^\ — (111 — O'ER + O-RL + o-rr)^ 
SCa = —C^\aLL ~ O-LR — O-RL + O-Rr)-, 
5C'^ = C^^\ — aLL + 0,LR — O.RL + 0,Rr) 

6Cs = gs{ALL + ALR + ARL + ARR), 

SC's = - 9si-ALL - Alr + Arl + Arr), 

6Ct = '2gT{aLL + orr), 

5C't = -2gT{-aLL + aRR). (19) 

The parameters dji, otji and Aji describe contributions to the low energy Hamiltonian 
from current-current interactions in terms of j-type of leptonic current and z-type of quark 
current. For example, cllr is the contribution to the Hamiltonian from left-handed leptonic 
current and right-handed quark current normalized by the size of the Standard Model (left- 
left current) interactions, gs and gT are formfactors at zero-momentum transfer in the 
nucleon matrix element of scalar and tensor currents. For more details, see the paper 
It should be noted, that 5Ci + 5C[ involve left-handed neutrinos and 5Ci — 5C[ is related to 
right-handed neutrino contributions in corresponding lepton currents. The analysis of the 
three experiments under consideration (a, A and B coefficient measurements) in terms of 
sensitivities {K~^) to aji, aji and Aji parameters is presented in the table IHl For the sake 
of easy comparison the sensitivities in this table are calculated under assumptions that all 
parameters [aji, aji and Aji) have exactly the same value, 1 ■ 10^^. The expected values of 
these parameters vary over a wide range from 0.07 to 10^® (see table ITTTl and, paper for 
the comprehensive analysis). The numerical results for the coefficients K in the table can be 
linearly re-scaled for the parameters (iij, aji and Aij in the range from 10^^ to 10~^ with an 
accuracy better than 10%. The relative statistical errors presented in the Table demonstrate 
discovery potentials of different experiments to new physics in terms of parameters ciij, aji 
and Aij. It should be noted, that the parameter cilr cannot provide sensitive information 
on new physics at the quark level, unless we obtain the axial-vector coupling constant gA 
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from another experiment, since in correlations diR appears in a product with gA (see j7[). 
For discussion of significance of each of these parameters to models beyond the standard 
one see 0. 

TABLE II: Relative statistical error (K) of the standard experiments to different types of inter- 
actions from new physics (ajj constants) provided that these constants have the same values of 
1 • 10-3. 





aLL 


aiR 


a-RL 


O-RR 


All 


Alr 


Arl 


Arr 


O^LL 


orr 


a 


0.17 


0.25 


135 


487 


1.43 


1.43 


283 


283 


0.19 


79 


A 


1.53 


0.63 


423 


1026 


13.1 


13.1 


860 


860 


1.82 


223 


B 


0.58 


1.21 


89 


347 


0.72 


0.72 


958 


958 


0.37 


59 



It should be noted the results in the tables H] and [H] are calculated with the estimated 
value of the parameter (a/(27r) Cy = 0.02. Numerical tests show that a change of this 
parameter by a factor two leads to changes of results in the tables by about 1%. 

The calculated integral sensitivities of different experiments to a particular parameter 
related to new physics can be used for the estimation of the experimental sensitivity when 
the experimental statistics is not good enough. For the optimization of experiments it is 
useful to know how manifestations of new physics contribute to the energy spectrum of 
the measurable parameter. As an example, the contributions from (Ilr, clrl and clrr to 
the spectra for the a, A and B correlations are shown on figures (0) - (jHl)- For uniform 
presentation all graphs on the figures are normalized by Nf = Cj^^lVudl"^ J f{E)dE, where 
f{E) is a{Ep), A{Ee) and A{Ef,), correspondingly. One can see that these contributions 



TABLE III: Possible manifestations of new physics 



Model 


L-R 


Exotic Fermion 


Leptoquark 


Contact interactions 


SUSY 


Higgs 


aRL 


0.067 


0.042 










aRR 
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0.01 








All + Alr 








0.01 


7.5 • 10-^ 


3 • 10-6 


Arr, + Arl 








0.1 






-All + Alr 






3 • 10-6 








Arr — Arl 






4 • 10-4 
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FIG. 1: Manifestation of ULR-type interactions on the a coefficient. 
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FIG. 2: Manifestation of URL-type interactions on the a coefficient. 
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FIG. 3: Manifestation of aLR-type interactions on the A coefficient. 



11 




FIG. 4: Manifestation of orl -type interactions on the A coefficient. 

<5A(aRR) 





FIG. 6: Manifestation of aiR-type interactions on the B coefficient. 



12 




FIG. 7: Manifestation of a^jL-type interactions on the B coefficient. 




FIG. 8: Manifestation of a_R_R-type interactions on the B coefficient. 

have different shapes and positions of maxima both for different model parameters and for 
different angular correlations. This gives the opportunity for fine tuning in the search for 
particular models beyond the Standard one in neutron decays. 

Using the approach developed here one can calculate the exact spectrum for a given 
model. For example, manifestations of the Left-Right model {aju = 0.067 and ajui = 0.075) 
in the measurements of the A and B coefficients are shown in solid lines on figures 1^ and ITUl 
The dashed lines show contributions from recoil effects and radiative corrections (without 
Coulomb corrections) assuming that (a/(27r) e(5) = 0.02. From these plots one can see 
the importance of the corrections at the level of the possible manifestations of new physics. 
The figure ^2 shows how these corrections for the coefficient B affected by the value of the 
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FIG. 9: Contributions from radiative and recoil corrections (dashed line) and from the left-right 
model (solid line) to the A coefficient. The curves are explained in the text. 
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FIG. 10: Contributions from radiative and recoil corrections (dashed line) and from the left-right 
model (solid line) to the B coefHcient. The curves are explained in the text. 

parameter {a/{27r) Cy) related to nuclear structure: dashed-doted, dashed and solid lines 
correspond to 0.01, 0.02 and 0.03 values for the parameter. 

We presented here results of analysis for only a number of parameters dij to illustrate a 
different level of sensitivities of experiments to the parameters. For the complete analysis 
of future experiments all aij, OLij and Aij parameters should be analyzed with a specific 
experimental conditions taken into account. 
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FIG. 11: Contributions from radiative and recoil corrections to the B coefficient for (a/(27r) Cy) = 
0.01 (dashed-doted line), (a/(27r) e^) = 0.02 (dashed line), and (a/(27r) e^) = 0.03 (solid line). 

IV. CONCLUSIONS 

The analysis presented here provides a general basis for comparison of different experi- 
ments of neutron /3-decay from the point of view of the discovery potential for new physics. 
It is also demonstrates that various parameters measured in experiments have quite differ- 
ent sensitivities to the detailed nature of the (supposed) new physics and can, in principle 
be used to differentiate between different extensions to the Standard Model. Thus neutron 
decay can be considered as a promising tool to search for new physics, which may not only 
detect the manifestations of new physics but also define the source of the possible deviations 
from predictions of the Standard model. Our results can be used for optimization of new 
high precision experiments to define important directions and to complement high energy 
experiments. Finally we emphasize that the usual parametrization of experiments in terms 
of the tree level coefficients a, A and B, is inadequate when experimental sensitivities are 
comparable or better to the size of the corrections to the tree level description. This is 
expected in the next generation of neutron decay experiments. Therefore, such analysis is 
needed for these experiments. One has to use the full expression for neutron beta-decay in 
terms of the coupling constants. In other words, the high precision experiments should focus 
on the parameters important for physics rather than on the coefficients a, A and B which 
are sufficient only for low-accuracy measurements. 




Ee (MeV) 
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